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^ ' Abstract. We studied the impact of a dynamical threshold on the /-/curve— the 

/^S ' relationship between the input and the firing rate of a neuron— in the presence of 

^^ ^ background synaptic inputs. First, we found that, while the leaky integrate-and- 

^vj , fire model cannot reproduce the /-/curve of a cortical neuron, the leaky integrate- 

and-fire model with dynamical threshold can reproduce it very well. Second, 
we found that the dynamical threshold modulates the onset and the asymptotic 
behavior of the /-/ curve. These results suggest that a cortical neuron has an 
r ) ^ adaptation mechanism and that the dynamical threshold has some significance 

>_y , for the computational properties of a neuron. 
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1. Introduction 

Neurons receive thousands of synaptic inputs, which are then transformed into output 
spike trains. A common assumption in neuroscience is that the firing rate— the average 
number of spikes per unit time— of a neuron conveys information and is used for 
computation in tlie brain. Tliis assumption is supported by a number of experimental 
studies [U [21 [16] . It is therefore very important to study the relationship between the 
input current and the firing rate (/-/ curve) in the presence of background synaptic 
inputs in order to understand the computation performed by a neuron. 
Recently, the leaky integrate-and-fire with dynamical threshold (LIFDT) model [HI 
[SJ [HI dH dH [131 [13 has been proposed so as to include spike frequency adaptation. 
Although the LIFDT model is simple, it can reproduce important properties of cortical 
neurons such as the negative interspike interval correlation O [15] and the /-/ curve 
of pyramidal neurons driven by a fluctuating current J13j . It is also reported that 
the spike threshold of a cortical neuron in vivo is not constant but depends on the 
preceding spike times [TU] . 

Our understanding of the effect of the dynamical threshold on the computational 
properties of a neuron is lacking. We studied whether the leaky integrate-and-fire 
(LIF) and LIFDT neuron can reproduce the /-/ curve of a cortical neuron [7] and 
how the dynamical threshold modulates the /-/ curve in the presence of background 
synaptic inputs. 

2. Method 

2.1. Leaky integrate-and-fire neuron 

Here, we briefly introduce the leaky integrate-and-fire (LIF) neuron [12]. The 
membrane potential V of the neuron obeys the first-order differential equation 

c'^ = ^9i.{V-V^)+I{t), (1) 

where C is the membrane capacitance, ^l is the leak conductance, Voo is the resting 
potential of the neuron, and I{t) is the input current. When the membrane potential 
reaches the threshold 6*0, a spike is generated and we instantaneously reset V{t) to the 
resetting potential V^. Wc adopted the parameters of the LIF neuron from Chance et 
al. [7], where C = Q.^[nF/cTn\ ,9l = 0.025[AiS/cm2], V^o = -65[mV], V, = -60[mV], 
and6lo = -54[mV]. 

2.2. Leaky integrate-and-fire neuron with dynamical threshold 

Wc briefly introduce here the LIFDT neuron [H H] [HI [TU [H [H [H] . The membrane 
potential V of the neuron obeys Eq. ([1]). In the following, we take the parameters 
C, (7l, Voa, Vy to be the same as in the LIF neuron. When the membrane potential 
reaches the threshold 9{t), a spike is generated and we instantaneously reset V{t) to 
the resetting potential Vi. The emission of a spike causes the threshold to increase by 
an amount Ae and then decay to its resting value 9oo exponentially: 

du a — Uqq , X — > „, , ,„, 

^^-^— +A.5:^(t-...), (2) 

k 

where Tg is the time constant of the dynamical threshold and t^ is the fc-th spike 
time. The sum is taken over all the spikes generated by the neuron up to time t. 
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The initial condition of the thresliold is its resting value: 6{Q) = Ooo- The parameters 
of the dynamical threshold are 9oo = — 54[mV], adopted from Chance et al. [T]., and 
Tg — 80 [ms], adopted from Liu and Wang [T?| . 

2.3. Synaptic input 

We considered an input current I{t) = m + Is{t), which is the sum of a constant 
current m and a background synaptic current /g. The synaptic current is described 

by 

Ut) =. gEit)iVE - Vit)) + fnit){Vi - Vit)), (3) 

ffE.i(i) =y]aE,iexp ^ , 

where ^e,! are the excitatory (E) and the inhibitory (I) synaptic conductances, Ve,i 
are their respective reversal potentials, te,! are their respective time constants, oej 
are their respective peak conductances, and i^.' are the fc-th spike times of the 
respective presynaptic neuron. The spike times of the respective presynaptic neuron 
are generated by an independent homogeneous Poisson process with the same rate 
7- The parameters are te = 5[ms], tj = 10[ms], Ve = 0[mV], Vi = — 80[mV], 
Oe = 0.01[yuS/cm^], and ai — 0.04[yLtS/cm^], adopted from Chance et al. [7]. 

2.4. f-I curve 

For the LIF and LIFDT neurons, we calculated the /-/ curves, which give the 
relationship between the constant current m and the firing rate given a fixed synaptic 
input rate 7. The firing rate / is the number of spikes emitted by a neuron per unit 
time, / = Nsp/T, where Ngp is the number of spikes and T is the observation time 
interval. To calculate the /-/curve, we simulated Eqs ([l]), 121), and ([3]) with time step 
St~ 0.01 [ms] and the time interval r= 50[s]. 

3. Results 

3.1. Irreproducibility of the f-I curve by the LIF neuron 

We compared the /-/ curves of a cortical neuron obtained by Chance et al. [7] (Figure 
lA) to those of the LIF neuron (Figure IB). We found that the LIF neuron cannot 
reproduce two main features of the /-/ curves of a cortical neuron. The first is the 
onset of the /-/ curve. The onset of the /-/ curves of a cortical neuron are linear, 
whereas those of the LIF neuron are nonlinear. The second irreproducible feature is 
the asymptotic behavior of the /-/ curve. The asymptotic behavior of the /-/ curves 
of a cortical neuron are sublinear, whereas those of the LIF neuron are linear. 



3.2. Effect of the dynamical threshold on the f-I curve 

We investigated the effect of the dynamical threshold on the onset of the /-/ curve. 
To quantify the nonlinearity of the onset, we fitted the onset of the /-/curve by using 
the power function 

h{x) = ci'x' + Co, 
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Figure 1. Comparison of tlie /-/ curves of a cortical neuron and the LIF neuron. 
(A) The /-/ curves of a cortical neuron with varying levels of synaptic input 
rate (adapted from Chance et al. [7] with kind permission from the author and 
Elsevier). (B) The /-/curves of the LIF neuron. Synaptic input rate 7 = 67.5[Hz] 
(closed diamonds), 7 = 135[Hz] (open circles), 7 = 270[Hz] (closed squares), 
7 = 405 [Hz] (open triangles). 



where [3 represents the nonhnearity of the onset, and cq, ci are the parameters. Figure 
[2j\ shows the onset of the /-/curves of the LIFDT neuron with three values of Ag while 
the synaptic input rate 7 is kept fixed. Figure [2j3 shows the dependence of the onset 
nonhnearity (3 on Ag. The onset of the LIF neuron [Ag = 0) is nonlinear—/? = 2.0— 
while the onset of the LIFDT neuron with large Ag is linear-like: /3 = 1.3. The 
dynamical threshold linearizes the onset of the /-/ curve. 

Next, we investigated the effect of the dynamical threshold on the asymptotic behavior 
of the /-/ curve. Here, we derived the asymptotic formulae of the /-/ curve. We used 
the diffusion approximation |131 I14j , wherein the excitatory and inhibitory synaptic 
conductances are approximated by the diffusion processes, and then we neglected the 
time correlation of the synaptic conductances. 



g^{t) « aETE7 + aE\/ — CeC^), .9i(i) 



T-El , 



aiT"i7 + fli 



^6(i) (4) 



where ^e, i{t) are independent Gaussian white-noise processes of zero mean and unit 
SD. Using Eqs ([I]), ©, and (g]), wc obtain 



C^ = -gl{V-V^)+m + s{V)m (5) 



where 



5l = .9l + 7(aETE + am), 

s\V) - ^-{alMV - V^f + a\n{V - Vif}, 

Here, ^(i) is a Gaussian white-noise process of zero mean and unit SD, and we assume 
that the synaptic input is balanced, (cete + aiTi)Voo = ceteI^e + clitiVi. We write 
the firing rate of the LIF neuron (Eq. ([5]) ) with the threshold 6'o as /lif("^, s{V), Oq), 
and the firing rate of the LIFDT neuron as /dt("t-, s(y)). For large constant input 
s{V) <^ m, we can neglect the contribution of the noise s{V), 

/lif(to, 3,6*0) « /lif(to, 0,6*0), /dt(to, s) « /DT(m,0). 
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The firing rate of tire LIF neuron can be written as 



9l 






CiOo - Vr) 



const. 



(6) 



The firing rate of tire LIFDT neuron is given by the self-consistent solution [13] of 

/DT(m, 0) ^ /lif( m, 0, Oo, + AgTgfMm, 0) ). (7) 

From Eq. ([7]), we can obtain 

Vid.oo -V,)^+AmAgTg/C- 



/dt(»7i, 0) 



-Vr] 



2AgTg V C^eTe 

Thus we can obtain the asymptotic formula of the /-/ curve for large m, 

m 



const. (8) 
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Figure [3|\ shows the /-/ curve of the LIF neuron, and Figure [3j3 shows the /-/ curves 
of the LIFDT neuron. In the large m, these /-/curves approach Eqs ([6]) and ([8]). The 
dynanrical threshold changes the asymptotic behavior of the /-/curves: the asymptotic 
behavior of the LIF neuron is linear (/ ~ m), whereas that of the LIFDT neuron is 
sublinear (/ ^ to^/^). 
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Figure 2. The effect of the dynamical tliresliold on tlie /-/ curve. (A) Tlie /-/ 
curves of the LIFDT neurons fitted by the power function (lines). Parameters: 
7 = 135[Hz] (fixed) while Ag = (circles), Ag = 0.5 (triangles), Ag = 3 (crosses). 
(B) The nonlinearity f3 of the /-/ curve is plotted as a function of Ag . Parameters: 
7 = 135 [Hz] (fixed). 



3.3. Reproducibility of the f-I curve by the LIFDT neuron 

We compared the /-/ curves of a cortical neuron obtained by Chance et al. [7j (Figure 
4A) to those of the LIFDT neuron (Figure 4B). The LIFDT neuron can reproduce 
the two main features of the /-/ curve of a cortical neuron: the linear onset and the 
sublinear asymptotic behavior. Our results are summarized in Table [TJ 
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Figure 3. The effect of the dynamical threshold on the asymptotic behavior 
of the /-/ curve. (A) The /-/ curve of the LIF neuron. The solid hne is the fit 
of Eq.®. Parameter: 7 = 135[Hz]. (B) The /-/curve of the LIFDT neurons. 
Solid lines are the fit of Eq.JSj. Parameters: 7 = 135[Hz] (fixed) while Ag = 1.0 
(circles), Ag = 5.0 (crosses). 
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Table 1. Main features of the /-/ curves 
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Figure 4. Comparison of the /-/ curves of a cortical neuron and a LIFDT 
neuron. (A) The /-/ curves of a cortical neuron with varying levels of synaptic 
input rate (adapted from Chance et al. [7] with kind permission from the author 
and Elsevier). (B) The /-/ curves of a LIFDT neuron. Parameters: Ag = 5.0 
(fixed), synaptic input rate 7 = 67.5[Hz] (closed diamonds), 7 = 135[Hz] (open 
circles), 7 = 270[Hz] (closed squares), 7 = 405[Hz] (open triangles). 
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4. Discussion 

4.I. Conclusion 

We studied the effect of a dynamical threshold on the /-/ curve. It is confirmed that 
the dynamical threshold is essential for reproducing the /-/ curve of a cortical neuron 
and that it modulates the onset and the asymptotic behavior of the /-/ curve. These 
results suggest that a cortical neuron has an adaptation mechanism and that this 
significantly infiuences the computation of a cortical neuron. 

4-2. Unresolved problems 

Here we state two unresolved problems relating to our study. The first is the 
clarification of the effect of adaptation on the network behavior of neurons. It 
would be interesting to study the effect of adaptation on not only the behavior 
of a neuron but also on that of interacting neurons. The second open problem is 
the clarification of the effect of the adaptation mechanism on the Hodgkin— Huxlcy- 
type model neurons [1]. For the sake of simplicity, we studied LIF and LIFDT 
neurons in this paper. However, Hodgkin— Huxley-type models are known to be more 
realistic neuron models. It is currently thought that the adaptation mechanism mainly 
arises from M-type currents, mAHP-type currents, and slow sodium currents in the 
Hodgkin— Huxley- type model [31 [51 [T^] . It would be interesting to study the effect of 
these currents on the /-/ curve. 
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